ON TAILS OF PERPETUITIES 



PAWEL HITCZENKOt 



Abstract. We establish an upper bound on the tails of a random variable that arises 
as a solution of a stochastic difference equation. In the non-negative case our bound is 
similar to a lower bound obtained by Goldie and Griibel in 1996. 



1. Introduction 
A random variable R satisfying the distributional identity 
(1) R = MR + Q, 

where (M, Q) are independent of R on the right-hand side and = denotes the equality in 
distribution, is referred to as perpetuity and plays an important role in applied probability. 
The main reason for this is that it appears as a limit in distribution of a sequence {Rn) 
given by 

Rn = MnRn-l + Qn, n> 1, 

provided that limit exists (here, (M„, Q„) is a sequence of i.i.d. random vectors distributed 
like (M, Q) and Ro could be an arbitrary random variable; for convenience we will set 
Rq = 0). Systematic study of properties of such sequences was initiated by Kesten in [S] 
and they continue till this day. Once the convergence in distribution of (-Rn) is established, 
at the center of the investigation is the tail behavior of R. There are two distinctly different 
cases: 

P(|M|>1)>0 and P(|M| < 1) = 1. 
The first results in R having a heavy tail distribution, that is 

P{\R\ > x) ~ Cx~\ 

for a suitably chosen constant n and some constant C (see the original paper of Kesten 
[5] or [2]), while in the second case the tails of R are no heavier than exponential. This 
was observed by Goldie and Griibel in |3]. Some subsequent work is in [3], but the full 
picture in this case is not complete. The purpose of this note is to shed some additional 
light on that case by establishing a universal upper bound on the tails of \R\. In a special, 
but important, situation when Q and M (and thus also R) are non-negative our bound is 
comparable to a lower bound obtained by Goldie and Griibel in [3]. 
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2. Bounds on the tails 

For a random variable M such that \M\ < 1 and < 5 < 1 define ps '■= P{1 — S < 
\M\ < 1). Then, as has been shown in [3] (see also the equation (2.2) in [4j) if < M < 1 
and Q = q {q being a positive constant), then for < c < 1 and x > q we have 

Since ln(l — cq/x) < —cq/x, for any particular value of c, say c = 1/2, this immediately 
gives 

N / ln(l - c) , , ,21n2 , 

P{R > X) > exp( x\n{pcg/x)) = exp( a; lnpg/(2x)). 

cq q 

Our aim here is to supply an upper bound of a similar form. While our result does not 

give the asymptotics of P{R > x) as x — )■ oo, it shows that it essentially behaves like 

exp{^x\npc2g/x) for some positive constants ci,C2. Specifically, we prove 

Proposition 1. Assume \Q\ < q, \M\ < 1, and let R be given by ([T]). Then, for sufficiently 
large x 

P{\R\ > x) < exp(— xlnp2g/x)- 
Thus, ifQ = q>OandO<M<l then 

/21n2 X / X /I 

exp(— — xlnpg/(2x)) < P{R > X) < exp{—x\np2g/x)- 

Proof. If -P(|M| = 1) > then, as was proved in [3j, R has tails bounded by those of an 
exponential variable, so we assume that \M\ has no atom at 1. Fix < 6 < 1 and define a 
sequence (T^) as follows 

To = 0, r^ = inf{A;> 1: |Mr_,+fc| < 1-5}, m > 1. 

Then T^'s are i.i.d. random variables, each having a geometric distibution with parameter 
1 — ps- Furthermore, \Mk\ < 1 — 5 if = Ti + ■ ■ ■ + Tj for some ^ > 1 and \Mk\ < 1 
otherwise. Therefore, 

m 

Y[\Mk\<{l-Sy for Ti + --- + T,- <m<Ti + --- + T,- + T,-+i. 
fc=i 

This in turn implies that 

fc-1 fc-1 

I E n ^^^j ^ E n i^^-i < ^1 + (1 - + (1 - 6fT, + . . . = ^(1 - 6f-'n. 

k>l j=l k>l j=l fc>l 

Therefore, ii \Q\ < g we get 

(2) P{\R\ >x)< PiY, n l^.l ^ -) ^ PiT.^k{l - S)'-' > -). 

fc>l j=l ^ k>l ^ 
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To bound the latter probability we use a widely known argument (our calculations follow 
[H proof of Proposition 2]). First, if T is a geometric variable with parameter 1 — p then 



e^{l — p) e 



X 



p l-^(eA_i 



j=i j=i ' i-p 

provided e^p < 1. Thus, writing t in place of a;/g in the right-hand side of ([2]), for A > 
we have 



P(E(1 - (5)^-'Tfc >t) = P(exp(A ^(1 - 5)^-'Tfc) > e^*) < e-^'Ee^^^>^^''^^-^^" . 

k>l k>l 

If A satisfies e'*'p < 1 then pe'^^^"^-''' ^ < 1 for every A; > 1 as well, and by independence of 
(Tk), the expectation on the right is 

k=l l-p^"^ k=l l-p^ > 

Now, choose A > so that ^{e^ - 1) < |. Then, as 1/(1 - u) < e^^ for < m < 1/2, for 
every /c > 1 we get 

1 - 



1 - 1^(6^(1-^)'=-^ - 1) - l-p 

Therefore, the rightmost product in ([3]) is bounded by 

p 



< exp(2^— (e-'^^-'^^ - 1)). 



exp(2-^5:(e^(-^)^-^)-l)). 
1 — V ^-^ 



'■ k>l 

We bound the sum in the exponent as follows 

fc>l i>l j>l •'' k>l 

_ ^ A-'' 1 1 ^ A-'' _ - 1 

j>i ^ ^ j>i J 

Combining the above estimates we get that 

(4) Pm > qt) < exp(-tA + ^ + — ). 

1 — p d 

provided that A satisfies the required conditions, that is: 

e^p < 1 and j^i^^ - 1) < ^• 

Clearly both are satisfied when e^p < 1/2. 

We finish the proof by making a suitable choice of A. Since we are assuming that \M\ 
has no atom at 1 and we are interested in large x, we may assume that S is small enough 
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SO that Ps < 1/3. This condition imphes that 2ps/{l — ps) < ^Ps so that the last term in 
the exponent of is bounded by 3ps{e^ — Now let t = 2/6. Then (jl]) becomes 

Pi\R\ > qt) < exp(-A| + - + J^£^) < exp(-i(A - 3p,(e^ - 1))). 

1 — Ps 

Set A = In(gi^) so that e^ps = |. This choice of A is within the constraints and maximizes 

the value of A — 3ps{e^ — 1), this maximal value being 

ln(-^) - 3ps{:^ - 1) = In(-) - (1 + ln3) + 3^^ > ]- Hl/ps), 
^Ps ^Ps Ps ^ 

with the inequality valid for sufficiently small ps (less than e~^/9 for example). Thus, using 
t = 2/6 we finally obtain 

P{\R\ > qt) < exp(-^ In(lM)) = exp(^ Inpa/t), 

or, in terms of x, 

X 

Pi\R\ > x) < exp(— lnp2g/x)- 

□ 
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